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Abstract 



Recent advances in the QCD corrections to b — ► s'y decay in the MSSM in- 
clude i.) evaluation of the relevant operators, Wilson coefficients and anoma- 
lous dimension matrix elements for the various MSSM effective theories valid 
at scales beyond Q = Mw, H-) calculations of most of the needed anomalous 
dimension matrix elements to next-to-leading order for scales J$ Q < Mw , 
and Hi.) calculations of 0(a s ) virtual and bremsstrahlung corrections to the 
b — > S7 decay operators at scale Q ~ mj. We assemble all these known 
results to gain an estimate of B{b — > S7) for the parameter space of the 
minimal supergravity model (mSUGRA). We find a much reduced scale de- 
pendence of our result compared to usual leading-log evaluations. Comparison 
with the latest CLEO results yields stringent constraints on parameter space. 
Much of mSUGRA parameter space is ruled out for fi < 0, especially for 
large tan (3. We compare these results with other constraints from cosmology 
and non-standard vacua. Also, we compare with expectations for discovering 
mSUGRA at LEP2, the Tevatron and the CERN LHC. 
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I. INTRODUCTION 



Particle physics models including weak-scale supersymmetry (SUSY) are amongst the 
most compelling candidates for physics beyond the Standard Model (SM). Of this class of 
models, the minimal supergravity model (mSUGRA) stands out as at least the most popular 
framework for performing searches for SUSY, and can justifiably be called the paradigm 
model for weak scale supersymmetry ||. The mSUGRA model can be characterized briefly 
by the following attributes 

• particle content and gauge symmetries of the Minimal Supersymmetric Standard 
Model (MSSM), i.e. a supersymmetrized version of the two-Higgs doublet SM, plus 
all allowed soft SUSY breaking terms, 

• a desert between the weak scale and the unification scale, which allows for gauge 
coupling unification, 

• universal boundary conditions for soft-SUSY breaking terms mo, mi/2, Aq and Bq 
implemented at the unification scale Mx, 

• electroweak symmetry is broken radiatively, and is a consequence of the large top 
quark Yukawa coupling. 

The various weak scale parameters are related to GUT scale parameters via renormalization 
group equations (RGE's). Typically, all weak-scale sparticle masses and mixings are then 
determined by the parameter set 

m o, m i/2, A , tan/5 and sign(fi), (1.1) 

along with the measured value of m t . We take m t = 170 GeV throughout this paper. 

Of course, not all values of the above parameter set are allowed. For some values, 
electroweak symmetry is not broken appropriately. For other values, a charged or colored 
sparticle may be the lightest SUSY particle (LSP), in conflict with cosmology and searches 
for exotic nuclei and atoms. In addition, there exist constraints from negative searches for 
sparticles at the Fermilab Tevatron pp and LEP2 e + e~ colliders. In particular, we note the 
recent bound that > 79 GeV for a gaugino-like lightest chargino ||. Additionally, 

in the absence of /^-violating interactions, there exist bounds on parameter space from the 
relic density of neutralinos produced in the Big Bang [BJ; for certain regions of parameter 
space, the neutralino relic density Qh 2 > 1, which implies a universe younger than 10 billion 
years, in contradiction at least with the ages of the oldest stars in globular clusters. Finally, 
several recent papers have mapped out regions of mSUGRA parameter space where there 
exist vacua deeper than the standard minimum |J. These non-standard vacua constraints 
may not be rigorous, however, if one entertains the notion that the universe may have settled 
into a false vacuum. 

The parameter space of the mSUGRA model (as well as many other models 0) may 
also be constrained by data from rare meson decays, such as the branching fraction for 
B — ► Xs'-f. A recent analysis by the CLEO collaboration finds for the inclusive decays 
B(B -> X s j) = (2.32 ± 0.67) x 10" 4 , with 95% CL upper and lower limits (including 
systematic errors) of 4.2 x 10~ 4 and 1 x 10 -4 , respectively |§. Such data, when compared 
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against theoretical predictions of b — > s"y, have been shown to be very restrictive for both 



two-Higgs doublet models (2HDM) || and supersymmetric models [10,11 • For a type II 
2HDM, loop contributions involving the top quark and charged Higgs add constructively 
with SM loops involving top and W bosons. In the MSSM, there exist other contributions 
involving squark-chargino loops, squark-neutralino loops, and squark-gluino loops. The 
latter two contributions are much smaller than squark-chargino loop contributions, and are 
frequently neglected. The squark-chargino loop contribution can add either constructively 
or destructively with the W and H + loops, leading to allowed or forbidden regions of SUSY 
model parameter space. 

These calculations are usually performed by evaluating lowest order matrix elements of 
effective theory operators at a scale Q ~ m&. All orders approximate QCD corrections are 
included via renormalization group resummation of leading logs (LL) which arise due to a 
disparity between the scale at which new physics enters the b — > sj loop corrections (usually 
taken to be Q ~ M w ), and the scale at which the b — > sj decay rate is evaluated (Q ~ ra^). 
The resummation is most easily performed within the framework of effective field theories. 
Above the scale Q = Mw (all scales Q ~ M w are equivalent in LL perturbation theory), 
calculations are performed within the full theory. Below Q = Mw, heavy particles are 
integrated out of the theory, leading to an effective Hamiltonian 

AC 8 

Heff = -4WE^(e)Oi(Q)> (1-2) 

V2 i=l 

where matching between the two theories occurs at Q = M w and yields the values of 
Ci(Q = Mw)- In Eq. O], the Ci(Q) are Wilson coefficients evaluated at scale Q, and the Oj 



are a complete set of operators, given, for example, in Ref. [12|. Resummation then occurs 



by solving the renormalization group equations (RGE's) for the Wilson coefficients 

Q^C i {Q)= lji C j {Q), (1.3) 
where 7 is the 8x8 anomalous dimension matrix (ADM), and 

7= £ 7p,+( £ )V " + — (L4) 

The matrix elements of the operators Oi are finally calculated at a scale Q ~ my, and 
multiplied by the appropriately evolved Wilson coefficients to gain the final decay amplitude. 

The LL QCD corrections just described yield enhancements in the b — > 57 decay rate of 
factors of 2 — 5 |13|,14]. The resulting LL calculation yields an answer which is ambiguous 



depending upon which precise scale choice is chosen for evaluation of matrix elements of the 
operators Oi. Variation of the scale ^ < Q < 2rrib yields approximately a 25% uncertainty 
in the theoretical calculation. This uncertainty provides the greatest source of error in 
currently available theoretical calculations [|14| | . To reduce the theoretical uncertainty, one 



must proceed to a next-to-leading-log calculation (NLL) of the b — > S7 decay rate. 

Recently, a number of theoretical developments have been made towards the goal of 
a NLL calculation of B(b — ► sj). Cho and Grinstein showed [T5| that if there are two 
significantly different masses contributing to the loop amplitude (such as m t and Mw), then 
in fact there can already exist significant corrections to the Wilson coefficients at scale Mw- 



3 



In this case, one must create an effective theory by first integrating out the heavy top quark, 
apply RGE running between the scales m t and M w , and then integrate out the W boson to 
arrive at the operator set in |1.2|. The above procedure gives a ~ 20% enhancement to the 



SM value of B(b — > 57). In supersymmetric models, many more heavy particles are present. 
Anlauf has shown how to perform a similar analysis for the case of the MSSM ||16| . These 
corrections are considered to be next-to-leading order effects. 

In addition, various terms of the ADM in Eq. [T3| have been calculated to next-to-leading 
order (NLO). The 0(a 2 s ) terms 7^- for i, j = 1 — 6 have been calculated by Ciuchini et. al. 



|T7| 1, while the corresponding terms for z,j = 7,8 are given by Misiak and Miinz Of 
the remaining terms which mix the i = 1 — 6 and j = 7, 8 Wilson coefficients, only 7I7 is 
relevant, but its evaluation requires a three loop calculation. A preliminary report on the 
calculation of 727* indicates that it is only a small effect ||19[ |. 



Finally, the QCD corrections to the operators Oi must be included [20,21]. Recently, 



Greub, Hurth and Wyler have reported results |2T[ on the complete virtual corrections to 
the relevant operators O2, and Og. Combining these with the bremsstrahlung corrections 
1 22|,21] results in cancellation of associated soft and collinear singularities. A combination of 
the QCD corrected operator matrix elements < s~f\Oi\b > with the complete 0(a s ) corrected 
Wilson coefficients at scale Q ~ nib will result in a NLL calculation of the B(b — > S7) decay 
rate. 

In this paper, we have nothing new to add to the calculational procedure for evaluating 
B(b — > 57). Our goal in this paper is to bring together the above pieces of a NLL calculation 
of B(b — > 57) and to interface with the mSUGRA model so that detailed comparisons to 
data can be made in parameter space. In so doing, we simply neglect the missing piece of 
the calculation, 727* . This will result in some small scheme dependence of our results, and 
in some additional scale dependence, so our calculation will not be truly NLL. Hence, we 
label it as QCD improved, in hope that 727 turns out small, as preliminary reports suggest. 

Our goal as well is to evaluate the B(b — > 57) rate as a function of mSUGRA parameters 
and compare with the recent CLEO results, to find favored or excluded regions of parameter 



space. We find the b — > S7 constraint to be really very strong, as noted previously [10,11 



We compare the b — > sj results to other recent results on relic density constraints || and 
regions of non-standard vacua ||. Finally, we note the effect of B{b — > sj) on expectations 
for discovering mSUGRA at various collider experiments. Toward these ends, in Sec. 2 
we present various details of our QCD improved calculation for B(b — > 57). In Sec. 3, we 
report on our results as a function of mSUGRA parameter space, and compare with other 
constraints and expectations for collider searches. We summarize in Sec. 4. 



II. CALCULATIONAL METHOD 

In this section, we outline our procedure for calculating B(b — > 57) as a function of 



mSUGRA parameter space. Our first step of course is to input the parameter set |1.1| and 
solve for the superparticle masses and mixings via running of the mSUGRA RGE equations 
between Mz and Mqut and imposing the appropriate minimization criteria using the 1-loop 
corrected effective potential. We iterate the running back and forth between the two scales 
six times using the usual Runge-Kutta method; this results in a convergent spectrum of 



superparticle masses. The procedure is described more fully in Ref. |[23|| . On the last run 
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down but one, we take note of the various superparticle masses. On the final run from Mqut 
to Mz, we implement the procedures outlined in Anlauf [16] to ultimately obtain the needed 
Wilson coefficients C 7 (M W ) and C 8 (M W ). In addition, the values of Ci(M w ) for % = 1 — 6 
are given in Ref. [|17] to 0(a s ). 

Given a heavy particle of mass mn and a light particle of mass contributing to a 
b — ► S7 loop, the procedure of Anlauf is as follows. First, at scale Q = mj, the heavy particle 
is decoupled from the theory and the corresponding effective field theory is constructed. The 
leading terms of the expansion of Q in terms of x = {^) 2 are calculated and evolved to 
scale Q = using the ADM including both QCD and electroweak interactions. Taking 
only leading terms in x restricts the operator basis to dimension 5 and 6 operators ||16|| . At 
Q = rriL, the remaining part of C*, which has been evolved down to using only the EW 
ADM, is added together with eventual contributions coming from decoupling of the lighter 
particle in the loop. As a last step, the equations of motion are applied to obtain Cj(Mw) 
and Cs(Mw)- In the case that ttll < Mw, the evolution is stopped at Mw- In practise, 
this is almost never a problem, in light of the new bound > 79 GeV for gaugino-like 

charginos from LEP2. In our calculation, we include contributions from tW, tH~ and Wi<jj 
loops, for i = 1 and 2, and q~j = u L , ur, Cl, Cr, t\ and t 2 , and neglect contributions from 
Zi<jj and gq~j loops, which should have much smaller contributions. 

As an illustration of these results, we show in Fig. 1 the final value of the Wilson 
coefficient Cj(Mw) from the various loop contributions just discussed, as a function of mi/ 2 , 
where mo = mu2, tan/? = 2 and Aq = 0. In frame a) \i < 0, while in frame 6), \l > 0. The 
tW SM contribution is just ~ —0.23, and of course doesn't vary versus SUSY soft-breaking 
parameters. In both frames, the tH~ contribution is negative, and decreases in absolute 
value as mi/2 increases, since the value of m#- is increasing. Note this contribution adds 
constructively to the SM tW contribution; when combined, these give the large constraints 
on type II 2HDM's |§. For /i < 0, most of the SUSY particle loop contributions are also 
negative in this case, which leads to the significant constraints to be given in Sec. 3. The 
exceptions are the large positive contributions from \V2t2 and W\q, which cancel some of 
the large negative contributions. Alternatively, for \x > 0, we see in frame b) that there are 
several positive as well as negative contributions to CV(Mvy). In this case, one can achieve 
rates for B(b — > sj) which are equal to or even smaller than the SM value. 

The next step in our calculation is to implement the NLO ADM to calculate the running 
of the Cj's from M w down to Q ~ m b . The terms 7}- for i, j — 1 — 6 have been calculated 



in Ref. [17]. The corresponding ADM elements for i,j = 7,8 have been calculated in Ref. 



T8| . The remaining NLO ADM elements for 7^- for i = 1 — 6 and j = 7, 8 have not yet 



been published. The most important of these is 7.^, since C^Mjy) — 1 while the remaining 
Ci(Myy) which mix into C-j are ~ pl| . Preliminary results of the three-loop calculation of 



727^ indicate it is small [TjJ, so in our calculation, we take 727 = 0. 

We show in Fig. 2 the evolution of the set of Wilson coefficients Cj for % = 1 — 8, from 
their values at Q — Mw down to Q — 1 GeV. Frame a) shows the evolution including just 
LO terms in the ADM, while frame b) includes the NLO ADM contributions mentioned 
above. We see that in general, the NLO effects are small. The exception is for Cr, for which 
the NLO correction differs from LO by ~ 10%. 

Finally, we must include the matrix elements < s^f\Oi(Q)\b > at NLO. We neglect the 
3 , 4 , 5 and Oq contributions, since the corresponding Cj are small, as shown in Fig. 
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2. Furthermore, the matrix element < s'yjOx |& > is exactly zero. The bremsstrahlung 



graphs have been calculated in Ref. [22[ for b — > s^g, while the complete virtual corrections 



to the operators O2, O7 and 0$ have been calculated in Ref. pi| . We implement these 
results, which ensure the proper cancellation of infrared and collinear singularities. Since 
our final results are not completely NLL, we will have some remaining scheme dependence. 
All our calculations have been performed within the naive dimensional regularization (NDR) 
scheme, in which the calculational building blocks have been given. We neglect throughout 
our calculation any long distance effects p4] on the b — > sj decay rate. 
Our final numerical result is given by 

B(b^sy) = ^B sh (2.1) 

1 si 

where 

r(6 -> 57) = T virt + T brem . (2.2) 



In the above, T virt is given by Eq. (5.6) of Ref. and T brem is given in Ref. [p2] , pl}| , while 



T s i is given by Eq. (5.9) of Ref. pifl . Numerically, we take the combination jy-^r = 0.95 
and B sl = 0.104. 

The inclusion of the various above mentioned QCD improvements leads to a result for 
B{b — > sj) which has significantly reduced scale-dependence. We illustrate the scale depen- 
dence of our result in Fig. 3, where we plot B{b — > sj) versus Q, where ^ < Q < Imy,. 
The LL calculation for the SM value is denoted by the dashed curve, which yields 
B{b — > sj) = 2.9 d= 0.7, or a 25% uncertainty due to scale choice. The QCD-improved 
result is shown by the solid curve. In this case, the prediction is B(b — > sj) = 3.2 ±0.3, and 
the error due to scale choice uncertainty is reduced to ~ 9%. The CLEO measured central 
value is denoted by the solid horizontal line. The la limits are denoted by dotted lines, and 
the 95% CL limits are denoted by solid lines. We see that the SM prediction lies somewhat 
above the CLEO measured result, although it is well within the 95% CL region. 



III. NUMERICAL RESULTS AND IMPLICATIONS 

A. Results for B(b — > S7) 

We present our main results on the b — > 57 branching ratio as contours of constant 
branching fraction in the mo vs. rrii/2 plane. This allows a direct comparison of previous 
work on mSUGRA constraints and also collider expectations to be made with the present 
work. In Fig. 4, we show the B{b — > sj) contours for A = for a) tan/3 = 2, [i < 0, b) 
tan/3 = 2, n > 0, c) tan/? = 10, /i < and d) tan/3 =10, /i > 0. Each contour must be 
multiplied by 10~ 4 . The values of B(b — > sj) shown are for Q = m b . The regions marked 
TH are excluded by theoretical constraints: either the lightest SUSY particle (LSP) is not 
the neutralino Zi, or electroweak symmetry is improperly broken. The regions marked EX 
are excluded by negative SUSY particle search experiments at Fermilab Tevatron or LEP. 
The new bound on > 79 GeV on gaugino-like charginos from LEP2 is indicated by the 
dashed contour. 
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In frame a), we find B(b — > sj) to be large for small values of mo and m^. This is due 
to constructive interference amongst many of the SUSY and SM loop contributions, as can 
be gleaned from Fig. la. As m and m 1( / 2 increase, the SUSY particles and charged Higgs 
boson all increase in mass, and their loop contributions become small: thus the value of 
B(b — > 57) gradually approaches its SM value in the upper-right region of each frame. We 
note the CLEO 95% CL upper bound on the inclusive rate for B(b — » 57) is 4.2 x 10 -4 , so it 
becomes evident that a significant region to the lower-left of frame a) will become excluded. 
In frame b), for fi > 0, the results are significantly different. In this case, there are many 
interfering loop contributions (see Fig. 16), so the B(b — ► 57) rate is much closer to the SM 
value, and can even drop below it. Since all the contours lie within the CLEO 95% excluded 
band, this frame remains unconstrained by B(b — > 57). For frame c), with tan/3 = 10 and 
/i < 0, we find very large values of B(b — > 57) throughout the entire region of the plane 
shown. Almost all of the plane shown gives values of B(b — > 57) greater than the CLEO 
95% CL bound, and so will be excluded. Values of mi/ 2 > 500 GeV are required to reach 
an allowed region for this choice of mSUGRA parameters. Finally, in frame d), we show 
the B(b — > S7) result for tan/3 = 10 and fi > 0. As in frame b), there exists substantial 
interference amongst the various loop contributions. The interference in this case is so great 
that a large fraction of the plane actually has B(b —>■ 57) values below the SM value. The 
B(b —>■ sj) rate increases with mi/ 2 to reach its SM value. 

Fig. 4 showed results for the GUT scale trilinear coupling Aq = 0. Changing the value of 
Aq will change the weak scale A parameters, which can result in changes to the top squark 
mass matrix and mixing angles. This can then affect the W{tj loop contributions. To show 
the effect of changing Aq, we show in Fig. 5 the B{b — > 57) contours in the mo vs. A 
plane, for fixed mi/2 — 200 GeV, and all other parameters as in Fig. 4. There are some 
small TH excluded regions in the left-hand corners of the tan/3 = 2 frames. In frame a), we 
see that the B(b —>■ 57) rate varies mainly with mo rather than Aq, and that only a small 
portion of the plane is above the CLEO 95% excluded value of B(b — > 57) = 4.2 x 10~ 4 . In 
frame b), the B(b — > 57) rate varies slowly versus parameters, and the entire plane shown 
is allowed. Frame c) again varies slowly with Aq, and is entirely exluded. Frame d) has 
significant variation against parameters, but is still entirely allowed. 

B. Relationship to Other Constraints on the mSUGRA Model 

Our next task is to compare the constraints from b — > 57 with other constraints, and 
derive conclusions relevant for collider searches. Towards this end, we show in Fig. 6 
regions of the m vs. mi/ 2 plane which are excluded by CLEO data on Bib — > sj) at 
95% CL. We match against the theoretical result from this paper. To obtain the excluded 
region, the B(b — > 57) rate must fall outside the CLEO allowed values for all scale choices 

< Q < 2m,ft. The relevant excluded region for frame a) lies to the lower- left of the solid 
contour labelled b —>■ 57. We show as well the dashed contour, below which m^. < 79 
GeV, in violation of recent LEP2 chargino searches The region to the left of the line 
of open circles is where non-standard minima of the mSUGRA model scalar potential lie 
0. This region may not be truly constrained if one is willing to accept that our universe 
may have settled into a false vacuum. The region to the right of the solid contour labelled 
by Qh 2 is excluded by cosmological considerations 0. In /^-conserving models where the 
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lightest neutralino is the LSP, it can be an excellent candidate to make up the bulk of the 
dark matter of the universe. Such a cold dark matter (CDM) particle would have been 
produced in abundance in the early universe, and their present day abundance can then be 
straightforwardly calculated. If the relic density is too high, then the calculated lifetime 
of the universe is too short to be consistent with various astrophysical observations. The 
contour denoted by Qh 2 denotes where the relic density Qh 2 = 1; for higher values of Qh 2 , 
the universe would be younger than 10 billion years old, and so is certainly excluded. In 
addition, a number of cosmological observations (COBE data, nuclear abundances, large 
scale structure) favor a universe formed with a 2:1 ratio of CDM to hot dark matter (HDM- 
e.g. neutrinos). This is called the mixed dark matter scenario (MDM). This cosmologically 
favored region lies between the dot-dashed contours, for which 0.15 < Qh 2 < 0.4. 

We combine the above mentioned constraints all on one plot, as in Fig. 6. We see that 
in frame a), if one combines the false vacua region with the B(b — > 37) constraint, much 
of the cosmologically favored MDM region is ruled out! Recent calculation of parameters 



associated with fine-tuning in the mSUGRA model |25j actually prefer the small 7% and 
mi/2 regions that are excluded by B[b — > 57) in this case. 

In frame b) of Fig. 6, the entire region shown is allowed by B(b — > 57), and in fact 
would be favored over the frame a) results by the CLEO central value, which lies somewhat 
below the SM B(b — ► 57) prediction. Meanwhile, much of the region shown in frame c) 
is excluded by b —>■ 57, including all of the cosmological MDM prefered region. Finally, 
in frame d), the entire region is allowed by B(b — > S7). In fact, in this case, the region 
around mi/2 ~ 200 GeV actually agrees with the central value of the CLEO measurement 
of B(b — > 57), and overlaps considerably with the cosmological MDM favored region! The 
corresponding excluded regions from the false vacuum constraint |J and Bib — > 57) for the 
m vs. A plane for mi/ 2 = 200 GeV are shown in Fig. 7. 



C. Implications for Collider Experiments 

Next, we wish to draw some conclusions for future searches for mSUGRA at colliding 
beam experiments. Expectations for mSUGRA in the same m vs. m\i% planes have been 



worked out for the CERN LEP2 e + e collider [26 , various Fermilab Tevatron pp collider 



options P7| , and for the CERN Large Hadron Collider (LHC) a pp collider expected to 
operate at 14 TeV. 

We return attention to Fig. 6 a. If the mSUGRA model represents nature with parameters 
as in Fig. 6 a, then the B(b — > 57) exclusion region wipes out most of the parameter space 
accessible to LEP2 SUSY searches. There are two exceptions |]26[| . There is a small region 
with mo ^ 200 GeV and m 1( / 2 — 100 GeV where charginos could be accessible to LEP2 
operating around ^fs ~ 190 GeV. In this region, the lightest neutralino Z\ can still be a 
good CDM candidate, since m^ ~ ^If-- thus, relic neutralinos can annihilate away via a Z 
boson pole in the s-channel. The other possibility for LEP2 is to discover a light Higgs boson 
h in the region beyond the B(b — > sj) exclusion contour. Much of the discovery reach of the 
Fermilab Main Injector (/ Cdt = 2 fb^ 1 ) upgrade is also wiped out by the b — > 57 constraint 
for the region m ^ 180 GeV. As with LEP2, there is some remaining region of accessibility 
for seeing clean trileptons from W\Z 2 — > 3£ around m ^ 200 GeV and mj/2 ~ 120 GeV. 
The TeV33 upgrade would be able to see in addition a large slice of parameter space beyond 
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the b — > 57 excluded region up to mi/ 2 ~ 275 GeV for m < 100 GeV, again via clean 
trileptons [^7|. LHC would of course be able to scan well beyond the entire plane shown, up 
to values of mi/ 2 ~ 700 — 1000 GeV with just 10 fb _1 of data. Since B{b — > 57) decreases 
with increasing mi/2 and m , however, the CLEO data prefer the region of parameter space 
accessible to LHC experiments, rather than the regions accessible to LEP2 and Tevatron, in 
contrast to preferences from fine-tuning p5f . 

Fig. 66 is entirely unconstrained by B(b — > sj). However, we do note that in this case 
the B(b — > 57) rate decreases with decreasing m 1( / 2 in the low mu2 region. Hence, the 
CLEO data actually prefer the regions accessible to LEP2 and Tevatron experiments, as do 
fine-tuning calculations. 

For Fig. 6 c, most of the plane shown is excluded by B(b — > sj), so if nature chose 
these parameters, then LEP2 and Tevatron upgrades would see nothing, and the discovery 
of SUSY would have to wait for LHC, which could access the very heavy sparticle spectra 
that lie at parameter space points beyond mi/ 2 ~ 350 GeV. 

The parameter space of frame 6c? is entirely allowed by B(b — > 37), and in fact the 
CLEO central value actually agrees with the parameter space region around m-1/2 ~ 200 
GeV. Since B{b — > 57) is decreasing below the SM value with decreasing m^, the low 
m.1/2 region accessible to LEP2 and Tevatron upgrades is dis-favored by data. The region 
around m ~ 125 and m 1 / 2 ~ 200 is favored by B{b — > 57) and by cosmology, and has no 
non-standard minima in the scalar potential. In this favored region, unfortunately, neither 
SUSY nor Higgs particles would be accessible to LEP2 experiments. Only a fraction of this 
region would be accessible to Tevatron upgrades via trilepton searches. However, the CERN 
LHC would enjoy huge supersymmetric signal rates in this region, and in addition, direct 
production of sleptons would be visible as well. 

IV. CONCLUSIONS 

We have performed a calculation of B(b — > 57) as a function of the parameter space 
of the mSUGRA model. In doing so, we have included several improvements over the 
usual leading- log treatment. We have included corrections to the Wilson coefficients due to 
multiple scales in both the SUSY and SM loop contributions. We have also included the 
published NLO terms in the ADM elements needed for running the Wilson coefficients from 
a scale Q = M w down to Q ~ m b . We have not included the correction to 727, so that our 
final result is not NLL. Preliminary results from the three loop calculation of 727'' indicate 
it is only a small effect. Finally, we have included the virtual and bremsstrahlung graphs 
in the evaluation of the < s-y\Oi\b > matrix elements. The combination of all the above 
elements leads to a B(b — > 57) calculation with reduced uncertainty due to choice of scale. 

We plot our results as a function of mSUGRA model parameter space, and compare 
against recent results from the CLEO experiment. The comparison leads to allowed and ex- 
cluded regions of mSUGRA parameter space. In particular, we note that for some mSUGRA 
parameter choices, the mSUGRA B(b —>■ 57) prediction agrees better than the SM one. The 
resulting constraints on parameter space are very strong, and indicate that large regions for 
fi < and for large tan j3 are excluded. We compare briefly with expectations for collider 
experiments at LEP2, Fermilab Tevatron, and CERN LHC. The B(b — > 57) constraint rules 
out much of parameter space that would have been accessible to LEP2 and Tevatron exper- 
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iments. However, we note that the region of parameter space most favored by B(b — > S7), 
cosmology and standard minima of the scalar potential, around rrio ~ 125 GeV, m 1( / 2 ~ 200 
GeV, tan f3 ~ 10 and /i > might be accessible to Fermilab MI or TeV33 searches for clean 
trileptons; if not, the discovery of SUSY would have to wait until LHC experiments are 
performed. 
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FIGURES 



FIG. 1. We plot the value of the Wilson coefficient CV(Mw) versus m^, where mo = m^, 
tan/3 = 2, t4o = 0. In a), we take \i < and in 6) we take /U > 0. The various loop contributions to 
Ct(Mw) are denoted on the plot. The label q refers to the sum over q = ul,ur,cl,cr contributions. 

FIG. 2. Evolution of Wilson coefficients d(Q) for i = 1 - 8 from Q = M w to Q = 1 GeV, 
for (mo, fn\j2i A), tan/3, sign(fi)= (100,100,0,10,-1), where masses are in GeV units. Frame a) 
shows evolution to LO, while frame b) shows NLO evolution (except 7^ = 0). 

FIG. 3. A plot of the SM branching ratio B(b — > 57) versus scale choice Q, where 
< Q < 2m&. We plot the LL result, and in addition, our QCD-improved result. The the- 
oretical error diminishes from ~ 25% to ~ 9%. We also plot the CLEO measured central value, as 
well as la and 95% CL limits on the experimental result. 

FIG. 4. Plot of contours of constant branching ratio B(b — ► 57) in the mo vs. m^ 2 plane, 
where ^o = and mt = 170 GeV. Each contour should be multiplied by 10~ 4 . Frame a) is 
for tan/3 = 2, fj, < 0, b) is for tan/3 = 2, /i > 0, c) is for tan/3 = 10, /x < and d) is for 
tan (3 = 10, /i > 0. The regions labelled by TH (EX) are excluded by theoretical (experimental) 
considerations. The dashed contour corresponds to the latest LEP2 limit of m^ > 79 GeV for a 
gaugino-like chargino. 

FIG. 5. Plot of contours of constant branching ratio B(b — > 57) in the mo vs. ^o plane, 
where m\/2 = 200 GeV and mt = 170 GeV. Each contour should be multiplied by 10 -4 . Frame 
o) is for tan/3 = 2, /j < 0, b) is for tan/3 = 2, /i > 0, c) is for tan/3 = 10, /i < and d) is for 
tan [3 = 10, /i > 0. The regions labelled by TH (EX) are excluded by theoretical (experimental) 
considerations. 

FIG. 6. Plot of contours of various constraints on the mSUGRA model in the mo vs. m 1( / 2 
plane, where ^o = and mt = 170 GeV. The frames are as in Fig. 4. To the left of the 
contour marked by open circles is the region where minima occur in the mSUGRA scalar potential 
that are deeper than the standard one. The region below the dashed contour is excluded by the 
LEP2 limit that m^ > 79 GeV. The region to the right of the solid contour labelled Qh 2 is 
excluded because the universe would be younger than 10 billion years old (Wi 2 > 1). The region 
between the dot-dashed contours is favored by the cosmological mixed dark matter scenario, where 
0.15 < flh 2 < 0.4. Finally, the region to the lower-left of the solid contour labelled b — > 57 in frame 
a) is excluded at 95% CL by the analysis of this paper. The entire region in frame b) is allowed 
by b — > 57, while almost the entire region in frame c) is excluded by b — > S7. Finally, the entire 
region shown in frame d) is again allowed by the b — > 57 constraint. 
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FIG. 7. Plot of contours of various constraints on the mSUGRA model in the mo vs. Aq plane, 
where mi/2 = 200 GeV and nit = 170 GeV. The frames are as in Fig. 5. To the left of the contour 
marked by open circles is the region where minima occur in the mSUGRA scalar potential that are 
deeper than the standard one. The entire regions in frames a), b) and d) are allowed by b — > sj, 
while the entire region in frame c) is excluded by b — ► s*y. 
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